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We present measurements of CP-violating asymmetries and branching fractions for the decays 
B"*" LJTT^ , — > ujK^ , and — *■ ujK". The data sample corresponds to 232 million BB pairs 
produced by e^e~ annihilation at the T(45') resonance. For the decay _B" ^ i-^Kg, we measure 
the time-dependent CP-violation parameters S = O.Sltg gg ± 0.02, and C = -0.55toll ± 0.03. 
We also measure the branching fractions, in units of 10~^, B{B^ —* lotv^) — 6.1 ± 0.7 ± 0.4, 
B{B+ -> ujK+) = 6.1 ± 0.6 ± 0.4, and luK°) = 6.2 ± 1.0 ± 0.4, and charge asymmetries 

Ach{B+ LO-K+) = -0.01 ± 0.10 ± 0.01 and Aah{B+ ujK+) = 0.05 ± 0.09 ± 0.01. 

PACS numbers: 13.25.Hw, 12.15.Hh, ll.30.Er 



Measurements of time-dependent CP asymmetries 
in meson decays through a Cabibbo-Kobayashi- 
Maskawa (CKM) favored b ccs amplitude 0, have 
firmly established that CP is not conserved in such 
decays. The efi'ect, arising from the interference be- 
tween mixing and decay involving the CP- violating phase 
P = arg i-VcdV;jVtdVtl) of the CKM mixing matrix i, 
manifests itself as an asymmetry in the time evolution of 
the B°B° pah. 

Decays to the charmless final states (j>K^, K'^ , 

T]'K°, t:°K°, /o(980)X°, and ujK° are aU b qqs pro- 
cesses dominated by a single penguin (loop) amplitude 
having the same weak phase /3 Q . CKM-suppressed am- 
plitudes and multiple particles in the loop complicate the 
situation by introducing other weak phases whose contri- 
butions are not negligible; see Refs. .5, 6] for early quan- 
titative work in addressing the size of these effects. We 
define AS as the difference between the time-dependent 
CP-violating parameter S (given in detail below) mea- 
sured in these decays and S — sin2/3 measured in char- 
monium K'^ decays. For the decay ujK^, these 

additional contributions are expected to give AS ~0.1 
0, Q , although this increase may be nullified when final- 
state interactions are included 8]. A value of AS" incon- 
sistent with this expectation could be an indication of 
new physics ,9] . 

We present an improved measurement of the time- 
dependent CP-violating asymmetry in the decay i?" — s- 
ujK'^, previously reported by the Belle Collaboration 
based on a sample of ~30 events [l^. We also mea- 
sure branching fractions for the decays — > ojK'^, 
£?+ cj7r+, and 5+ ojK^ (charge-conjugate decay 
modes are implied throughout), and for i?"*" — > o^tt^, 
and B^ —> uK'^ , we measure the time- integrated charge 
asymmetry Ach = (r^ -r+)/(r" +r+), where is the 
width for these charged decay modes. In the Standard 
Model Ach is expected to be consistent with zero within 
our experimental uncertainty; a non-zero value would in- 
dicate direct CP violation in this channel. 

The data were collected with the BABAR detector 
at the PEP-II asymmetric e+e" collider. An integrated 
luminosity of 211 fb^^, corresponding to 232 million BB 
pairs, was recorded at the T{AS) resonance (center-of- 
mass energy ^/s = 10.58 GeV). Charged particles are 
detected and their momenta measured by the combina- 
tion of a silicon vertex tracker (SVT), consisting of five 



layers of double-sided detectors, and a 40-layer central 
drift chamber, both operating in a 1.5 T axial magnetic 
field. Charged-particle identification (FID) is provided 
by the energy loss in the tracking devices and by the mea- 
sured Cherenkov angle from an internally reflecting ring- 
imaging Cherenkov detector (DIRC) covering the central 
region. A K/n separation of better than four standard 
deviations (a) is achieved for momenta below 3 GeV/c, 
decreasing to 2.5cr at the highest momenta in the B de- 
cay final states. Photons and electrons are detected by a 
CsI(Tl) electromagnetic calorimeter. 

From a B'^B^ pair produced in an T{AS) decay, we re- 
construct one of the B mesons in the final state / = i^K^, 
a CP eigenstate with eigenvalue —1. For the time evolu- 
tion measurement, we also identify (tag) the flavor (i?° 
or B'^) and reconstruct the decay vertex of the other B. 
The asymmetric beam configuration in the laboratory 
frame provides a boost of (3^ = 0.56 to the T(4S'), which 
allows the determination of the proper decay time differ- 
ence At = tf — ttag from the vertex separation of the two 
B meson candidates. Ignoring the At resolution (about 
0.5 ps), the distribution of At is 

e-|At|/r 

P(At) = \ItAw± (1) 

At 

(1 - 2w) {Ssm{AmdAt) - C cos{AmdAt))]. 

The upper (lower) sign denotes a decay accompanied by 
a P° (-6°) tag, T is the mean lifetime, Anid is the 
mixing frequency, and the mistag parameters w and Aw 
are the average and difference, respectively, of the proba- 
bilities that a true 5° {B'^) meson is tagged as a B'^ (-B"). 
The parameter C measures direct CP violation. If C = 0, 
then S = sin2/3 + AS". 

The flavor-tagging algorithm ^ has seven mutually ex- 
clusive tagging categories of differing purities (including 
one for untagged events that we retain for yield deter- 
minations). The measured analyzing power, defined as 
efficiency times (1 — 2w)^ summed over all categories, is 
(30.5 ± 0.6)%, as determined from a large sample of B- 
decays to fully reconstructed flavor eigenstates (Bflav)- 

We reconstruct a B meson candidate by combining a 
7r+, or with an a; Tr+Tr^Tr*^. We select K'^ 
TT+TT^ decays by requiring the 7r~''7r~ invariant mass to be 
within 12 MeV of the nominal K'^ mass and by requiring 
a flight length greater than three times its error. We 
require the primary charged track to have a minimum 
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of six Cherenkov photons in the DIRC. We require the 
7r+7r~7r° invariant mass (ma^) to be between 735 and 825 
MeV. Distributions from the data and from Monte Carlo 
(MC) simulations guide the choice of these selection 
criteria. We retain regions adequate to characterize the 
background as well as the signal for those quantities taken 
subsequently as observables for fitting. We also use in the 
fit the angle 9h, defined, in the w rest frame, as the angle 
of the direction of the boost from the B rest frame with 
respect to the normal to the ld decay plane. The quantity 
Ti = I cos Oh \ is approximately flat for background and 
distributed as cos^ Oh for signal. 

A B meson candidate is characterized kine- 
matically by the energy-substituted mass mES = 

\J (5S + Po • Pb)^/Eq — pg and the energy difference 

AE = E"^ — ^y/s, where (i?o,po) and {Eb,Pb) are 
four-momenta of the T{AS) and the B candidate, 
respectively, and the asterisk denotes the T{AS) rest 
frame. We require, assuming the — > wtt"*" hypothesis, 
\AE\ < 0.2 GeV and 5.25 < toes < 5.29 GeV. 

To reject the dominant background from continuum 
e'^e~ qq events {q = u,d,s,c)^ we use the angle Ot 
between the thrust axis of the B candidate and that of 
the rest of the tracks and neutral clusters in the event, 
calculated in the T{AS) rest frame. The distribution of 
cos Ot is sharply peaked near ±1 for jet-like qq pairs and 
is nearly uniform for the isotropic B decays; we require 
I cos6't| < 0.9 (0.8 for the charged B decays). 

From MC simulations of B'^B'^ and B^B^ events, we 
find evidence for a small (0.5%) BB background contri- 
bution for the charged B decays, so we have added a BB 
component to the fit described below for those channels. 

We use an unbinned, multivariate maximum-likelihood 
fit to extract signal yields and CP-violation parameters. 
We use the discriminating variables toes; AE, tostt, Ti., 
and a Fisher discriminant The Fisher discrimi- 

nant combines five variables: the polar angles with re- 
spect to the beam axis in the T(45) frame of the B can- 
didate momentum and of the B thrust axis; the tagging 
category; and the zeroth and second angular moments of 
the energy flow, excluding the B candidate, about the 
B thrust axis We also use At for the B" luK'^ 

decay, while for the charged B decays we use the PID 
variables and Tr-, defined as the number of standard 
deviations between the measured DIRC Cherenkov angle 
and that expected for pions and kaons, respectively. 

For the i3° decay we define the probability 

density function (PDF) for each event i, hypothesis j 
(signal and qq background), and tagging category c 
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V, (mEs^)P, iAE^)V, {r)V, (m^ J (2) 



extended likelihood function as 



(3) 



where Yj is the fit yield of events of species j, /j,c is the 
fraction of events of species j for each category c, and 
is the number of events of category c in the sample. We 
fix /sig,c to /bq^^.c, the values measured with the large 
Bflav sample The same likelihood function is used 
for the charged decays except that the hypothesis j also 
includes BB background, the tagging category is not 
used and the PDF is slightly different, involving flavor k 
(primary 7r+ or K^): 

r;, ^r,imEslv,iAEiTi)v,{r)r,{nii,)v,in^) . (4) 

The PDF VsigiAt, aAt,c), is the convolution of 
F{At; c) (Eq. ^ with the signal resolution function (a 
sum of three Gaussians) determined from the Sflav sam- 
ple. The other PDF forms are: the sum of two Gaus- 
sians for all signal shapes except Ti, and the peaking 
component of the tos^ background; the sum of three 
Gaussians for Vqq{At;c); an asymmetric Gaussian with 
different widths below and above the peak for Vj{J-') 
(a small "tail" Gaussian is added for Pqq{^)); Cheby- 
shev functions of second to fourth order for TC signal 
and the slowly- varying shapes of AE, m^T^, and TL back- 
grounds; and, for Pjiq (toes), a phase-space-motivated 
empirical function |l4j , with a small Gaussian added for 



BB 



(toes). 



where is the error on At for event i. We write the 



We determine the PDF parameters from simulation for 
the signal and BB background components. We study 
large control samples of i? ^ Dir decays of similar topol- 
ogy to verify the simulated resolutions in AE and toes, 
adjusting the PDFs to account for any differences found. 
For the qq background we use (toes, AE) sideband data 
to obtain initial PDF-parameter values but ultimately 
leave them free to vary in the flnal flt. 

We compute the branching fractions and charge asym- 
metry from fits performed without At or flavor tagging. 
The free fit parameters are the following: the signal 
and qq background yields (the BB yield, if present, is 
fixed); the three shape parameters of Vqq{J-); the slope 
of Vqq{AE) and Vqq{m3Tr)', the fraction of the peaking 
component of 'Pqq{m3Tr); ^ ^3; and, for the charged B 
decays, the signal and background Ach- 

Table HI lists the quantities used to determine the 
branching fraction. Equal production rates of B^B^ 
and B^B^ pairs have been assumed. Small yield biases 
are present in the fit, due primarily to unmodeled cor- 
relations among the signal PDF parameters. In Table 
in we include estimates of these biases, evaluated by fit- 
ting simulated qq experiments drawn from the PDF into 
which we have embedded the expected number of signal 
and BB background events randomly extracted from the 
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TABLE L Fit sample size, signal yield, estimated yield bias 
(all in events), estimated purity, detection efficiency, daughter 
branching fraction product, statistical significance including 
systematic errors, measured branching fraction, and corrected 
signal charge asymmetry. 



Quantity 




ujK+ 




Events in fit 


44175 


9145 


Signal yield 


274 ± 28 


266 ± 24 


100 ± 15 


Yield bias 


18 


16 


8 


Purity (%) 


34 


46 


46 


Eff. (e,%) 


21.8 


21.2 


23.0 




0.891 


0.891 


0.307 


exYlB. {%) 


18.2 


17.7 


6.4 


Significance (a) 


10.8 


13.0 


8.6 




6.1 ±0.7 


6.1 ±0.6 


6.2 ± 1.0 


Signal Ach 


-0.01 ±0.10 


0.05 ± 0.09 





fully simulated MC samples. The estimated purity in 
Table n is given by the ratio of the signal yield to the 
effective background plus signal, the latter being defined 
as the square of the error on the yield. 
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FIG. 1: The B candidate rriES and AE projections for B — > 
UJTV+ (a, h), B+ ^ LoK+ (c, d), and B° ujK° (e, f) shown 
for a signal-enhanced subset of the data. Points with error 
bars represent the data, the solid line the fit function, and 
the dashed line the background components. Note that the 
ljK'^ signal in the AE plot is displaced from zero since AE 
is defined for the ujtv^ hypothesis. 

Fig.nshows projections onto toes and AE for a subset 
of the data (including 45-65% of signal events) for which 
the signal likelihood (computed without the variable plot- 
ted) exceeds a threshold that optimizes the sensitivity. 

For the time-dependent analysis, we require \At\ < 20 
ps and a At < 2.5 ps. The free parameters in the fit are 
the same as for the branching fraction fit plus S, C, the 
fraction of background events in each tagging category, 
and the six primary parameters describing the At back- 




CP tag ' 



FIG. 2: Projections onto At for luK°. Data (points 

with errors), the fit function (solid line), background compo- 
nent (dashed line), and signal component (dotted line), for 
events in which the tag meson is (a) and (b) B", and (c) 
the asymmetry (A^^o — Ngo)/{Nga ± Ngo)- 

ground shape. The parameters t and Am^ are fixed to 
world-average values . Here we find a slightly smaller 
yield of 95±14 events and S = O.SltjJ j^, C = -O.bBto fe- 
The errors have been scaled by ~1.10 to account for a 
slight underestimate of the fit errors predicted by our 
simulations when the signal sample size is small. Fig. [21 
shows the At projections and asymmetry of the time- 
dependent fit with events selected as for Fig. ^ 

The major systematic uncertainties affecting the 
branching fraction measurements include the reconstruc- 
tion efficiency (0.8% per charged track, 1.5% per photon, 
and 2.1% per K*^) estimated from auxiliary studies. We 
take one-half of the measured yield bias (3-4%) as a sys- 
tematic error. The uncertainty due to the signal PDF 
description is estimated to be < 1% in studies where the 
signal PDF parameters are varied within their estimated 
errors. The uncertainty due to BB background is also 
estimated to be 1% by variation of the fixed BB yield 
by its estimated uncertainty. The Ach bias is estimated 
to be —0.005 ± 0.010 from studies of signal MC, control 
samples, and calculation of the asymmetry due to parti- 
cles interacting in the detector. We correct for this bias 
and assign a systematic uncertainty of 0.01 for Ach for 
both B+ UJTT+ and B+ ujK+. 



For the time-dependent measurements, we estimate 
systematic uncertainties in S and C due to BB back- 
ground and PDF shape variation (0.01 each), modeling 
of the signal At distribution (0.02), and interference be- 
tween the CKM-suppressed b — > ucd amplitude and the 
favored b cud amplitude for some tag-side B decays 
Hi (0.02 for C, neghgible for 5*). We also find that 
the uncertainty due to SVT alignment and position and 
size of the beam spot are negligible. The Sfiav sample is 
used to determine the errors associated with the signal 
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PDF parameters: At resolutions, taggin g effi ciencies, and 
mistag rates; published measurements |l5j are used for 
tb and Am^. Summing all systematic errors in quadra- 
ture, we obtain 0.02 for S and 0.03 for C. 

In conclusion, we have measured the branching frac- 
tions and time-integrated charge asymmetry for the de- 
cays [x)7r+ and ujK~^ and the branching 
fraction for B" ^ ujK°. Wc find B{B+ W7r+) = (6.1 ± 
0.7±0.4)xl0-6, B{B+ ujK+) = (6.1±0.6±0.4)xl0-6, 
B{B^ ujK^) = (6.2 ± 1.0 ± 0.4) x 10-^ Ach{B+ 
W7r+) = -0.01 ± 0.10 ± 0.01, and Ach{B+ ujK+) = 
0.05 ± 0.09 ± 0.01, where the first errors are statistical 
and the second systematic. These results are substan- 
tially more precise than earlier measurements jlTI] and a 
significant improvement over our previous measurements 
[18j, which they supersede. We also measure the time- 
dependent asymmetry parameters for the decay B*^ — > 
ujK°, S = 0.5lt°J^ ± 0.02 and C = -O.^bto fe ± 0.03, 
with a precision nearly a factor of two better than the pre- 
vious Belle Collaboration results 0|. If we fix C = 0, we 
find S = 0.60l^|. This value of S and the world-average 
value of sin2/3 yield a value of AS* = 0.12 ± 0.40, in 
good agreement with the expected value near zero. 
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